CREATED USING THE RSC LaTeX PCCP ARTICLE TEMPLATE - SEE www.rsc.org/electronicfiles FOR DETAILS 

ARTICLE TYPE www.rsc.org/xxxxxx | XXXXXXXX 



Localization and size distribution of a polymer knot confined in a 
channel 

Chihiro H. Nakajima, at Takahiro Sakaue/** 



o 

=3 



Received Xth XXXXXXXXXX 20XX, Accepted Xth XXXXXXXXX 20XX 
First published on the web Xth XXXXXXXXXX 200X 
DOI: 10.1039/b000000x 



We have examined the behaviors of a knotted linear polymer in narrow tubes using Langevin dynamics simulation to investigate 
the knot localization property in one-dimensional (ID) geometry. We have found that the knot is strongly localized in such a 
geometry. By observing the distribution function of the size of localized knot, we found the scaling behavior of the fluctuation 
around the most probable size with radius of confinement. Based on the analysis of the probability distribution of the knot size, 
q we show that the strong localization behavior and the fluctuation around the most probable size can be encompassed by a simple 
ry} argument based on the virtual tubes composed of parallel strands and overlapping among them. 
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1 Introduction 

Thanks to the advance in nanotechnology and single-molecule 
experiments, there has been a growing interest in the proper- 
ties of biopolymers in confined spaces. For instance, it is now 
possible to bring genomic-length DNA molecules into 100 nm 
size channels and to observe its static and dynamic behaviors 1 . 
For now, the existence of two different regimes depending on 
the channel size is well confirmed 2 , and such a basic knowl- 
edge is expected to provide a guide for the controlling and 
manipulation of single biopolymers in fabricated devices. 

In the present paper, we look into statistical properties of 
intramolecular entanglements formed in a confined polymer 
in narrow channels. Such topological defects are sometimes 
formed spontaneously in the course of experiments 3 , under- 
standing of which would be of practical as well as fundamen- 
tal importance. In addition, they can be made by design using 
manipulation techniques 4 , too. We shall mainly address the 
question about the size and fluctuation of such defects. 
\ The above question is directly related to the (de)localization 
phenomena of knots in closed ring polymers^. It is known 
that the phenomenon of knot localization is strongly related 
to the geometry of the space in which a polymer can explore 
its large scale conformation. Several recent studies have in- 
dicated that knots in three-dimensional (3D) space are likely 
to be weakly localized- - —. In 2D space (realized by the con- 
finement between parallel plates), the degree of the localiza- 
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tion is enhanced, and the strong localization is expected for 
flat knot s 11 ! 12 . On the other hand, polymers assume a com- 
pact globular state either in poor solvent conditions or in small 
closed cavities, where knots are known to be delocalized^. 



From this point of view, it is interesting to ask the behaviors 
of knots in ID geometries, i.e., a chain confined in a narrow 
tube. So far the study in this direction is not abundant. In par- 
ticular, there seems to be no investigation on the knot localiza- 
tion phenomena in ID channels, except for a related work by 
Mobius et al.— They studied the unknotting dynamics of an 
initially knotted linear semiflexible polymer in ID one chan- 
nels, and predicted the pathway, in which the macroscopic in- 
flation of knot takes place. Here and in the followings, we use 
the term "knot" for the defect (intramolecular entanglement) 
even for open linear polymers, unless confusion is expected. 

To clarify the static property of the knot in ID confine- 
ment, we have numerically investigated the behavior of a lin- 
ear chain containing a trefoil knot in a narrow tubeii. Al- 
though this trefoil may be eventually unknotted in the open 
chain with ends, one can naturally think of the precedent stage 
during which the knot exists in the identifiable manner. Our 
results indicate that knots are strongly localized in ID space 
and that the physical mechanism for the localization is entirely 
different from that in higher dimensions. In addition, the anal- 
ysis of the knot size distribution function reveals an important 
feature on the strand overlapping, which may be relevant to 
the statistics of ring polymers and star polymers confined in 
narrow tubes. 
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2 Models and simulations 




Fig. 1 (Color online)Setup of the simulation, (a): A typical 
conformation which appeares in the simulation with D = 6 and 
N = 500. (b): Illustration for the closure scheme. The gray dashed 
line corresponds to the supplemental arc. The part with lighter 
(orange) color is detected as the knotted region. 

We have performed Langevin dynamics simulations to in- 
vestigate the behavior of knotted polymers confined in a nar- 
row tube. We employ a linear chain model which consists of 
a sequence of beads of diameter a connected by springs: 



sp 



2a 2 L { 



(1) 



where ry+i is the distance between neighboring beads, the 
large modulus k sp /kBT = 60 (UbT being the thermal energy) 
keeps the bond length almost constant ~ a, and a is taken as a 
unity. A pair of beads interact by the repulsive Lennard- Jones 
potential: 

U int = (fcWfe) 12 - 1 } ,(2) 

I o ( r u > a ) 

where £ = 0.8&#7\ The polymer is confined in the cylindri- 
cal tube of radius D. By taking the tube axis as z axis in the 
Cartesian coordinate, the wall potential is implemented as, 

Vw = {^{te) 12 " 1 } (^*) ,0) 
I ( nj > a ) 



where r± t = yxf+yf is the distance of z-th bead from the 
tube axis. 

The motion of a model polymer chain in thermal and vis- 
cous medium is described by the Langevin equation; 



dt 2 



Y ltt 



dU 



-R/(0, 
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where m, y are mass and friction constant of monomers, 
U = U sp + Uint + Uw is the potnetial, and Rj(Y) is a Gaussian 
white noise on the z-th monomer, which satisfies the following 
properies, 

(R,-(0) = 0, (Rj(f)Ri(O) = 6yk B TSijS(t-t f ). (5) 

As an initial configuration, a trefoil loop is introduced by 
hand near the middle of chain in an arbitrary way. Starting 
from such a designed configuration, we equilibrate the system 
by numerically integrating the Langevin equation. Configu- 
rations were sampled after sufficiently long relaxation for the 
statistical analysis. Identification of the knotted component 
is based on calculations of the special value of the Alexander 
polynomial; A(a) with a = — 1. Since the Alexander polyno- 
mial is defined for a closed ring, we need a closure scheme 
to apply it to our linear chain. We label each monomer as 
n = 1 , 2, • • • , TV from one end, and make a closed ring by bridg- 
ing fti-th and ^2-th monomers with a supplemental arc (Fig. 
[Hb) ). Initially, both ends are bridged (n\ = 1, ri2 = N); then, 
the calculation of the Alexander polynomial for the knot dia- 
gram, which is obtained from the configuration of closed tra- 
jectory projected onto XZ plane, yields A(— 1) = 3, indicating 
that the trefoil knot is present in the range between n\ and 
ri2. We then increase (decrease) n\ (ni) one by one by cutting 
monomers from both ends, which is followed by the calcu- 
lation of the Alexander polynomial at each step. This itera- 
tive process is repeated until the point n\ = n\ and ri2 = n2, 
at which the value of the Alexander polynomial changes to 
A(— 1) = 1. This indicates that the closed trajectory obtained 
by bridging n\ and ti2 has a trivial topology. Then the knotted 
size is detected as A/& = ti2 — n\ . In other word, Nk is defined 
as the minimal component number yielding | A(— 1) | = 3 upon 
trial closing of the chain.The corresponding lateral size in ax- 
ial direction characterizes the spatial size of the knot (Fig. 
[TJa)). We note that the closure scheme generally has a finite 
probability to alter the topology at the stage of the closure us- 
ing an supplemental arc, which results in the misjudging of 
the size and the type of the knot—. Fortunately, such an error 
rate can be made very low in our quasi- ID system by utilizing 
its geometrical property. 

3 Results 

As described above, the simulation was started from a pre- 
pared initial configuration with a trefoil knot in the middle of 
the chain. When the inserted knot size Nk is large, we have 
observed that the system exhibits a spontaneous relaxation by 
expelling the strands toward the outside of the knotted part. 
It then settles in a stationary state in which the knotted part 
diffuses along the chain keeping its average size. A typical 
conformation in this stationary state is shown in Fig. Q](a) in 
the case ofD = 6,N = 500. 
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Fig. 2 (Color online)Dependence of (Lk) and (Nk) for various tube 
sizes on the total number of the monomers N. (Inset)Dependence of 
(L k ), (N k ) and N£ on D for the chain length N = 500. 



In Fig. [2 the averages of the knot size (Nk) and (Lk) 
are, respectively, plotted as functions of the total number of 
monomers N for various tube size D. Both Nk and Lk appear 
to be independent of the number of monomers of the polymer 
N, indicating the strong localization in ID geometry. 

To investigate the behaviors of confined knotted polymers 
in more detail, we plot in Fig. [3] the probability distribution 
function P(Nk) of the knot size in the stationary state for var- 
ious tube sizes. Here P(Nk) is defined and obtained as fol- 
lows. We first build the histogram H(Nk,Lk) where the pair 
values (NkXk) m ea ch sample is located Nk < Nk < Nk + AA^, 
and Lk < Lk < Lk + AL&. Here we adopt ANk = 1, AL& = 0.1. 
Then we approximate a two-dimensional density distribution 
of probability P(Nk,Lk) as 



P(Nk,L k ) : 



H(N k ,L k ) 



(6) 



(AAT^(AI*)(I> 44 ^,Z4))- 
From this density of probability, we construct the marginal 




Fig. 3 (Color online)Probability distributions of the knot size Nk 
with various tube sizes. 



distribution of Lk at fixed Nk as 



P(L k \N k 



P(N k ,L k ) 



f Lk P(Nk,L k )dL k 
ALkL N >MN' k ,L' k ) 



(7) 



and detect its peak position Z|(A^) for each value of N& We 
then calculate P(Nk) as P(Nk,L£). Because of the asymmetric 
profile of P(Nk), the peak value N% of the knot size is slightly 
smaller than the average value (Nk), but both exhibit qualita- 
tively the same dependence on D (see the inset of Fig. [2]). In 
Fig. HJa), we show typical snapshots of the knotted polymer 
confined in tube with D = 2 (left) and D = 4 (right), where the 
size of the knotted part is small (Nk < N£\ top), in average 
(N k ~ iV£; center), large (N k > N%; bottom). 

To quantitatively examine the spatial shape of the knotted 
part, we calculate the relative asphericity a and the relative 
prolateness %\ 




— A-2 — A3) (2A2 — A3 — X\) (2A3 — Ai — A2) 



2(Af- 



- A^ 



-A3 2 - 



Ai A2 + A2A3 + A3A1 



v3/2 



(9) 

where Ai > A2 > A3 are the eigenvalues of the gyration ten- 
sor of the knot. As shown in Fig. IHb), the ensemble of the 
knot shapes falls into a single curve in a a — 71 plane. The 
shape asymptotically approaches the rod limit (a — >► 1 , n — » 1) 
for large A^. The decrease in Nk from the large Nk limit is 
accompanied by the decrease in both a and %. This trend 
continues up to the most probable point Nk — N£, where 



This journal is ©The Royal Society of Chemistry [year] 



Journal Name, 20 1 0, [vol], 1 -0 | 3 




N/N k -1 ' ' N/N k -1 

Fig. 4 (Color online) (a)Typical confirmations of the knotted chain in the tube of size D = 2 (left) and D = 4 (right). Here shown are only the 
middle part of the whole chain containing the knotted part. (b)Parametric plot of the relative asphericity a and the relative prolateness K. 
(c)Dependence of a as a function of the normalized knot size Nk/N^ — 1. (d)Dependence of % as a function of the normalized knot size 
Njc/Nfc — 1. Error bars in (c) and (d) indicates the standard deviation for rough measure of fluctuations. 



(a, it) ~ (0.19,0.) <^Ai :A 2 : A 3 ~2: 1.5: 1. The value of a 
takes its minimum at N% . Further shrinkage of knot N < N£ 
leads to the slight increase in a, while % remains almost con- 
stant 7T = 0. While in very narrow tubes (D = 2, and 3), both 
a and n take slightly larger values at a given normalized knot 
size Nk/Nl — 1 [Fig. |3c) and (d)], these seem to converge 
to asymptotic curves for larger tube sizes. Note that in the 
region N <N£, each value of the prolateness has large fluctu- 
ation and Nk dependence of 7t has large ambiguity. But these 
dependence seem to have a tendency to approach n = as D 
becomes large. 

4 Discussion 

To understand the underlying mechanism of the strong knot 
localization and the behavior of the distribution function in 



Fig. [3j we consider the contribution of the free energy for 
confining the knotted polymer into the tube. The polymer can 
be divided into a knotted part with Nk monomers in the middle 
and two "arms" consisting of the N — Nk remanent monomers 
at both sides. Correspondingly we decompose the free energy 
F into F = Fk + F&. The arms can be pictured as a linear 
sequence of blobs of size D ~ ag v A , where v ~ 0.59 is the Flory 
exponent and g is the number of monomers per blob. The loss 
of degrees of freedom upon the confinement can be evaluated 
by assigning ~ ksT per blob 16 . We thus have 



F A 



N — Nk 

gA 



/ a \ !/ v 



(10) 



On the other hand, monomers in the knotted part feel more 
crowded environment. If we look at a perpendicular cross- 
section (X — Y plane) and count the number p of "primitive 
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Fig. 5 (Color online) (top)The schematic explanation of the concept 
of the virtual tube. (bottom)Functional form of the crowding factor 
y(D) as a function of D with A = 150. and B ~ 2.9. 




Fig. 6 (Color online)Collapse of the probability distributions with 
the scaling variables s = N k /N£-l and (P(N k ) /P(N£ ) J , which 
corresponds to q = 1 in Eq. Q3] The symbols are the same as those 
in Fig. [3] The inset represents the failure of the collapse with q = 0. 



path" strands [see Fig. 0top)], we find pa = 1 in arms and 
Pk > Pa in the knotted part. In the present case of the trefoil 
knot in a polymer, we have pt = 3. The section of polymer 
in the knotted part is thus effectively confined to a narrower 
tube. To implement such a feature, the simplest is to assume 
that each of p^ strand in the knotted part is confined in the 
virtual tube whose cross-sectional area is Sk = S/pk, where 
S~D 2 represents the cross-sectional area of the original tube. 
Then, the number of monomer per smaller blob due to the 
confinement is determined by the relation 



D 



ia 8k 



1/2 

with the crowding factor 7 = p k ' > 1 . This leads to the con 
finement free energy of the knotted part as follows; 



(11) 



W g k 



(12) 



Using Eq. (fTOb and (fT2li . the total free energy of the confine- 
ment can be represented as 



F 



k B T 



(^W'-D (13) 



where Fq ~ kBTN(a/D) l / y represents the conventional free 
energy of confining a linear polymer into tube, and the ad- 
ditional term is identified as the excess contribution due to 
the existence of the knot. This argument is valid for rela- 
tively large size of knot, and indeed explains the exponential 
tail in the knot size probability distribution (Fig. [3] and ) 
P(Nk) ~ ex P [—F(Nk)/kBT]. To reduce the crowded region as 
much as possible, the knot shrinks in ID confined geometry. 
The lower bound of Nk for this shrinkage is provided by the 
tube size, below which the strands in the knotted part should 
be more highly confined with its topological constraint 17 . We 
thus conclude that the knot size is determined by the tube size 
(Lk) c± D. Furthermore, if the optimum number of blobs in the 
knotted part n/gk has only weak or negligible dependence on 
the tube size £>, we expect the relation Nk — (Nk) ~ (D/a) l / v . 
These are consistent with the simulation results (Fig. [3]inset). 

The virtual tube assumption seems to succeed to explain 
the qualitative trend including the strong knot localization and 
the exponential tail in the knot size probability distribution. 
We now examine it on a more quantitative basis by look- 
ing at the crowding factor 7. To this end, we attempt to do 
a data collapse of the knot size distribution functions P(Nk) 
for various tube sizes by seeking for the scaling function for 
AF = F(N k ) - F(N£). According to Eq. ([13]), it is written as 



AF(s) 
k B T 



B(r l/v -l)s 



(14) 
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where s = (Nk/N%) — 1 is the normalized deviation from the 
optimum size, and a numerical coefficient B is introduced to 
express the relation between the optimum knot size and the 
tube size as N% = B(D/d) 1 / v . Replotting the numerical results 
in the form of Eq. (fT4l) , it turned out that the naive assumption 

1/2 

Y = p k = const, fails to collapse the data. Rather, we have 
found that the generalized form 



-As 



(15) 



with an exponent q ~ 1 provides a very good data collapse 
(See Fig. [6]). The value of A is roughly estimated as A ~ 
120. ~ 150. from the slope of the exponential tails of Fig. [6] 
Comparing Eqs. (ITU) and (fT5l) , we find that the crowding fac- 
tor should depend on the tube size as 



7(D): 



1 



/fl\?/v" 



(16) 



As shown in the bottom of Fig. [3 the crowding factor 7 de- 
creases with the increase in the tube size. This result could 
be interpretted as follows; while the most outer part of the as- 
sumed virtual tubes are composed of the original rather rigid 
wall potential, their inner parts originate from softer fluctuat- 
ing potential field due to the presence of other primitive path 
strands. This allows the partial overlap of the virtual tubes, 
the degree of which may increase with the tube size. The nar- 
rower the tube size is, the more the fluctuation of the primitive 
strands is suppressed. The D dependence of N£ shown in the 
inset of Fig. [2] determines B ~ 2.9. Using these independently 

determined value A and B in Eq. [161 Do ~ 4. is estimated as 

1/2 

the value at which 7(A)) = Pk w * m Pk = 3. The tube size Do 
signifies the very narrow space, inside which pk strands tightly 
fit to fill it. In fact, the slight deviations of the data points with 
D — 2 and 3 from the master curve in Fig. [4]and[6]are probably 
attributed to it. Although based on a rather rough estimation, 
such an internal consistency is regarded as a support for the 
line of our argument. 

5 Summary and perspectives 

There have been several recent studies on the (de)localization 
property of knots formed in a ring polymer. It has been ar- 
gued that, for flat knots (2D), the entropy of the looping is a 
relevant factor lending to the strong knot localization, the es- 
sential feature of which can be understood using the slip-link 
model 1 1 . In the present paper, we have shown that the knot is 
strongly localized in tubes (ID), too, but its physical mecha- 
nism is entirely different from that in 2D case. As might be 
evident from the analysis of the slip-link model, what is rel- 
evant here is not the entropy of the looping but the crowding 
due to the presence of parallel strands in the narrow tube. 



A useful way of taking such a crowding effect into account 
is to introduce the virtual tube concept, in which a plural num- 
ber pk of tubes are supposed to be present in accordance with 
the number pk of parallel strands. This naturally leads to the 
exponential probability distribution for the larger side of knot 
size, which has been well confirmed by numerical simulations. 

Our closer analysis of the knot size probability distribution 
has revealed that the crowding factor generally depends not 
only on the number pk of the parallel strands but also on the 
size of confining tubes. In addition to the knot localization 
phenomena, we expect that the present result would have an 
important consequence on the behaviors of circular polymer 
(pk = 2) which would be relevant to the bacterial chromo- 
some segregation— as well as the migration dynamics in nano- 
channels 20 . Another related system is a star polymer with / 
arms confined in tube, where we may set pk — fufr with //, 
f r = f — fi being the number of arms in left and right sides of 
the center 21 . 

Our main conclusion that the knot is strongly localized in 
ID geometry is in apparent contradiction with the prediction 
by Mobius et al. 13 . While the influence of the chain semi- 
flexibility included in the ref [12] still remains to be clarified, 
we expect that the presence of the crowding factor would be 
robust in ID geometries. This leads us to propose that the pre- 
dominant mode of the unknotting dynamics in ID geometries 
is not a knot inflation, but a diffusion of the localized knot 
along the chain. 
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